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Abstract. We show that the problem of deciding whether a polygonal knot 
in a closed three-dimensional manifold bounds a surface of genus at most g, 
is NP-complete. We also show that the problem of deciding whether a curve 
in a PL manifold bounds a surface of area less than a given constant C is 
NP-hard. 



1. Introduction 

In this paper we investigate the computational complexity of some problems in 
three-dimensional topology and geometry. We show that the problem of determin- 
ing a bound on the genus of a knot in a 3-manifold, is NP-complete. Using similar 
ideas, we show that deciding whether a curve in a metrized PL 3-manifold bounds 
a surface of area less than a given constant G is NP-hard. 

Determining whether a given knot is trivial or not is one of the historically cen- 
tral questions in topology. The problem of finding an algorithm to determine knot 
triviality was posed by Dehn Dehn's investigations into this area led to the 
formulation of the word and isomorphism problems, which played an important 
role in the development of the theory of algorithms. The first algorithm for the 
unknotting problem was given by Haken Haken's procedure is based on nor- 
mal surface theory, a method of representing surfaces introduced by Kneser |14j . 
Analysis of the computational complexity of this algorithm is more recent. Hass, 
Lagarias and Pippenger showed that Haken's unknotting algorithm runs in time at 
most c* , where the knot K is embedded in the 1-skeleton of a triangulated manifold 
M with t tetrahedra, and c is a constant independent of Af or if It was also 
shown in [Sj that the unknotting problem is in NP. 

The notion of genus was defined by Seifert in 1935 for knots in the 3-sphere, 
and extends directly to knots in an arbitrary 3-manifold M . Given a knot if, 
consider the class §(if ) of all orientable spanning surfaces for K. These are surfaces 
embedded in M with a single boundary component that coincides with K. Seifert 
showed that this class is non-empty for any knot K in the 3-sphere. For knots in 
a general manifold, §(if ) is non-empty when K represents a trivial element in the 
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first integer homology group of M. The genus g{K) of a knot K is the minimum 
genus of a surface in §{K), or oo if §{K) = 0. The genus measures one aspect of 
the degree of "knottedness" of a curve. 

The unknotting problem is a special case of the more general problem of de- 
termining the genus of a knot in a 3-manifold. Given a knot K in an orientable 
3-manifold and a positive integer g, this problem asks for a procedure to determine 
whether the knot genus of K , the minimal genus of an orientable spanning surface 
for X in a 3-dimensional manifold, is at most g. A knot is trivial, or unknotted, 
precisely when its genus is zero. We will show that the problem of determining 
the genus of a knot in a 3-manifold is NP-complcte. Previous results on this prob- 
lem were given in 0, where it was shown to lie in PSPACE, roughly the class of 
problems that run in polynomial space. No lower bounds on the running time were 
previously known. 

We work with 3-manifolds that are triangulated and orientable, and with ori- 
entable embedded surfaces. This is not a significant restriction, since all compact 
3-dimensional manifolds admit unique PL structures . A knot in a triangulated 
3-manifold M is a connected simple (non self-intersecting) closed curve in the 1- 
skeleton of M. Any smooth knot in a smooth manifold, or more generally any tame 
knot, is equivalent to a knot that lies in the 1-skeleton of some triangulation. 

We formulate the problem of computing the genus as a language-recognition 
problem in the usual way, see In 1961 Schubert in an extension of Haken's 
work, showed the decidability of the problem: 
Problem: 3-MANIFOLD KNOT GENUS 

INSTANCE: A triangulated 3-dimensional manifold M, a knot K in the 
1-skeleton of A/, and a natural number g. 
QUESTION: Does the knot K have g{K) < g? 

The size of an instance is measured by the sum of the number of tetrahedra in 
M. In Section Owe establish 

Theorem 1. 3-MANIFOLD KNOT GENUS is NP-hard. 

It was established in that 3-MANIFOLD KNOT GENUS is in PSPACE. 
We improve this bound in Section . 

Theorem 2. 3-MANIFOLD KNOT GENUS is NP. 

In combination these two results give: 

Theorem 3. 3-MANIFOLD KNOT GENUS is -complete. 

Theorem^is proved through a connection to ONE-IN-THREE SAT, a known 
NP-complete problem that will be reviewed in Section 3. The theorem carries 
out a construction that transforms an instance of ONE-IN-THREE SAT to an 
instance of 3-MANIFOLD KNOT GENUS. By "transform" we mean that an 
instance of one problem is changed to an instance of the second by a procedure 
that requires time polynomial in the size of the instance. To a boolean expression 
representing an instance of ONE-IN-THREE SAT we associate a positive integer 
g and a certain knot in a triangulated, compact 3-manifold. This knot bounds a 
surface of genus at most g exactly when there is a truth assignment to the boolean 
expression satisfying the requirements of ONE-IN-THREE SAT. Since ONE-IN- 
THREE SAT is NP-hard, this establishes that 3-MANIFOLD KNOT GENUS 
is also NP-hard. 
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In Section we prove Theorem El giving a certificate which demonstrates in 
polynomial time that a genus g knot K bounds a surface of genus at most g. The 
argument in 5 established that the unknotting problem is NP using the existence 
of a normal disk that lies along an extremal ray in the space of normal solutions, 
called a vertex surface in Jaco-Tollefson 11 . The existence of such an extremal 
normal surface of minimal genus spanning a knot is not known, so a new technique 
is needed (See 10 for known results here). This is provided in Theorem II 21 which 
gives an algorithm to count the number of orbits of a type of pseudogroup action on 
a set. Theorem 1121 seems likely to have more general applicability. In Section 0] we 
describe this algorithm and in Section |5l we apply it to a pseudogroup action that 
arises in the theory of normal surfaces. This allows us to determine in polynomial 
time the number of components in a normal surface described by an integer vector 
in Z^*. In particular we are able to certify that a normal surface is connected, 
orientable and has connected boundary. Since calculating the Euler characteristic of 
a normal surface can be done efficiently, establishing orientability and connectedness 
are the key steps in constructing a certificate of its genus. 

In Section El we extend the orbit counting algorithm to allow the counting of 
additional integer weight sums associated to each orbit. This allows for the polyno- 
mial time calculation of the genus of all the components of a fundamental normal 
surface, as well as a count of the number of components. 

The genus and the area of a surface are closely connected. In Scction[71we extend 
the methods developed in studying genus to study the problem of determining the 
smallest area of a spanning surface for a curve in a 3-manifold. We show that 
computing an upper bound on the area of a smallest area spanning surface is NP- 
hard. 

We refer to ^21 for a discussion of complexity classes such as NP and PSPACE 
and [25' for a discussion of complexity problems in low-dimensional topology. 
Remarks: 

(1) Knots are often studied in R'^ or rather than in a general manifold. 
Our methods show that determining knot genus in or S^, or any fixed 
manifold, is NP. It is not clear whether the corresponding problem remains 
NP-hard if one restricts consideration to knots in M.^ or S'^. 

(2) Casson has shown that a procedure to determine whether a 3-manifold is 
homeomorphic to the 3-sphere, following the 3-sphere recognition algorithm 
described in ^Hl and [22 j runs in time less than 3*p{t), where p{t) is a 
polynomial. In the direction of lower bounds, it was shown in that 
determining certain values of the Jones polynomial of alternating links is 
#P-hard. 

We are grateful to the referee for numerous suggestions on the exposition of this 
paper. 

2. Normal Surfaces 

General surfaces in 3-manifolds can wind and twist around the manifold in com- 
plicated ways. Kneser described a procedure in which surfaces can be "pulled taut" , 
until they take a simple and rigid position jl4| . In this normal position, they have 
very succinct algebraic descriptions. We use an approach to normal surfaces in 
triangulated 3-manifolds based on work of Jaco-Rubinstein |12j and Jaco-Tollefson 
[TT] . A normal surface 5 in a triangulated compact 3-manifold M is a PL-surface 
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whose intersection with each tetrahedron in M consists of a finite number of disjoint 
elementary disks. These are properly embedded disks that are isotopic to either 
triangles or quadrilaterals as shown in Figure ^ by an isotopy preserving each face 
of the tetrahedron. 




Figure 1. Elementary disks in a normal surface. 



Within each tetrahedron of AI there are four possible triangles and three possible 
quadrilaterals, up to a normal isotopy of M, an isotopy which leaves each cell of 
the triangulation of M invariant. W. Haken observed that a normal surface is 
determined up to such isotopies by the number of pieces of each of the seven kinds 
of elementary disks that occur in each tetrahedron, or a vector in ZJ^. A normal 
surface 5" is described by a non-negative integer vector v — v(5) G Z*"*, that gives 
the normal coordinates of S. There is a homogeneous rational cone Cm in R^* , called 
the Haken normal cone, that contains the vectors v(S') for all normal surfaces S in 
M. 

If V = {vi, W2, . . . , wyt) G MJ*^ , then the Haken normal cone is specified by linear 
equations and inequalities of the form 

Vii + — Vi^ + Vi^ (up to 6t equations) , 

v^>0 , for 1 < i < 7i . 
The first set of equations expresses matching conditions, which say that the number 
of edges on a common triangular face of two adjacent tetrahedra, coming from a 
collection of elementary disks in each of the tetrahedra, must match. For each trian- 
gular face there are three types of edges (specified by a pair of edges on the triangle) , 
which yield 3 matching conditions per face. Triangular faces in the boundary dM 
give no matching equations. The second set are called the positivity conditions. The 
cone Ca/ is rational, because the above equations have integer coefficients. We let 
Cm(Z) — Ca/HZ''* denote the set of integral vectors in the cone Cm- An additional 
set of conditions, the quadrilateral conditions, is required for an integral vector in 
the cone to correspond to the normal coordinates of an embedded surface. This 
condition states that of the three types of quadrilateral found in each tetrahedron, 
only one can occur in the vector with non-zero coefficient. The quadrilateral condi- 
tions are required because two distinct types of quadrilateral in a single tetrahedron 
necessarily intersect, and we are interested in embedded, non-self-intersecting sur- 
faces. A vector in the Haken normal cone that satisfies the quadrilateral conditions 
corresponds to an embedded normal surface. This surface is unique up to a normal 
isotopy. However it is important to note that the surface corresponding to a given 
vector in the cone Cm may not be connected. A normal surface meets each edge of 
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a triangulation in a finite number of points. Tfie sum of tliese intersection numbers 
over all the edges of a triangulation is called the weight of the normal surface. 
A fundamental normal surface is a normal surface S such that 

v{S) 7^ vi + V2 , with vi,V2 e eA/(Z)\{0} . 

In the terminology of integer programming, such a vector v(5') is an element of 
the minimal Hilhert basis M{Cm) of Cm, see Schrijver (50] Theorem 16.4]. A fun- 
damental normal surface is always connected, but connected normal surfaces need 
not be fundamental. A vertex minimal solution is a special kind of fundamental 
surface, one that corresponds to a solution of the normal surface equations that lies 
along an extremal ray of the cone of solutions and is not a multiple of another such 
extremal solution. Hass, Lagarias and Pippenger W, Lemma 6.1] gave a bound for 
the size of the vectors corresponding to any fundamental surface. 

Theorem 4. Let M be a triangulated compact 3-manifold, possibly with boundary, 
that contains t tetrahedra. 

• Any vertex minimal solution v 6 Z''* of the Haken normal cone Cm in M 
has max (vi) < 2^*^^ . 

• Any minimal Hilbert basis element v G Z*"* of the Haken fundamental cone 
Cm has max (vi) < t ■ 2'''*+^ . 

l<i<7t 

Schubert [21] showed that a surface of smallest genus spanning K can be found 
among the fundamental surfaces. 

Theorem 5. There is a minimal genus spanning surface for K which is a funda- 
mental normal surface. 

Similar to the theory of normal surfaces, though somewhat easier, is the theory 
of normal curves. These are curves on a surface that intersect each triangle in a 
collection of normal arcs, arcs that have endpoints on distinct edges. Normal curves 
arise as the boundaries of normal surfaces in a manifold with boundary. Since there 
are three such arcs in each triangle, normal isotopy classes of normal curves in a 
surface that contains t triangles are described by integer vectors in Z'^*. 

3. 3-MANIFOLD KNOT GENUS is NP-hard 

In this section we show how to reduce an instance of ONE-IN-THREE SAT 
to an instance of 3-MANIFOLD KNOT GENUS. Since ONE-IN-THREE SAT 
is known to be NP-hard, this establishes that 3-MANIFOLD KNOT GENUS 
is also NP-hard. The problem ONE-IN-THREE SAT concerns logical expres- 
sions involving collections of literals (boolean variables or their negations) gathered 
in clauses consisting of three literals connected with V's. The logical expression 
contains a collection of clauses connected with A's. 

Problem: ONE-IN-THREE SAT 

INSTANCE: A set U of variables and a collection C of clauses over U 
such that each clause c £ C contains 3 literals. 

QUESTION: Is there a truth assignment for U such that each clause in C 
has exactly one true literal? 

Schaefer (THI established that ONE-IN-THREE SAT is NP-complete. To prove 
Theorem □ establishing that 3-MANIFOLD KNOT GENUS is NP-hard, we 



6 



IAN AGOL, JOEL HASS, AND WILLIAM THURSTON 



show that an arbitrary problem in ONE-IN-THREE SAT can be reduced in 
polynomial time to a problem in 3-MANIFOLD KNOT GENUS. See Garey and 
Johnson [2] for a discussion and many examples of such reductions. 

Let U = {ui, U2, ■ ■ ■ Un} be a set of variables and C = {ci, C2, . . . Cm} be a set 
of clauses in an arbitrary instance of ONE-IN-THREE SAT. We will describe a 
knot K in a compact 3-dimensional manifold (with no boundary) and an integer 
g such that K bounds a surface of genus smaller or equal to g if and only if C is 
satisfiable so that each clause in C contains exactly one true literal. We construct 
the 3-manifold M in stages. First we construct a 2-dimensional simplicial complex, 
then we thicken this complex, replacing triangles with subdivided triangular prisms, 
getting a triangulated, 3-dimensional manifold with boundary, as indicated in Fig- 
ure 121 Finally we use a doubling construction, taking two copies of the manifold 
with boundary and gluing their boundaries together, to obtain a closed 3-manifold. 




Figure 2. A branching surface B, shown in cross-section, is 
"thickened" to produce a triangulated 3-manifold with boundary. 

To begin, we form a type of singular surface B, that we call a branching surface, 
by identifying boundary curves of a collection of 2n -I- 1 surfaces with boundary, 
each forming what we refer to as a piece of the branching surface. We construct this 
collection of surfaces as follows. Let ki be the number of times that the variable Ui 
appears in the collection of clauses, and ki be the number of times that the negation 
Ui of Ui appears. For i 1, . . . , n, let and Fu. be genus one surfaces with h + l 
and ki + 1 boundary curves respectively. Also set Fq to be a planar surface with 
n + m + 1 boundary curves. One of these boundary curves will later become the 
knot K. The branching surface B is constructed by identifying these surfaces along 
appropriate boundary components as indicated in Figure [3| 

Branching occurs when more than two boundary curves are identified along a 
single curve. We identify pairs of boundary curves by giving a homeomorphism 
between them. Up to isotopy, this is determined by specifying an orientation on 
the curves and setting the homeomorphism to be orientation reversing. We first 
fix an orientation on each of Fq, F^^ and Fn.. This induces an orientation on each 
boundary curve. All identifications will involve gluing a boundary component of 
Fui or Fu. to a boundary component of -Fbj and we require this gluing to be an 
orientation reversing homeomorphism. Label hy K,ui, . . . Un, ci, . . . €,„ the 1+n+m 
boundary components of Fq. The boundary component K of Fq, which will become 
our knot, has nothing identified to it. For each i, 1 < i < n, one boundary curve 
from the surface F^ is identified to Ui. The remaining ki boundary components of 
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Figure 3. A branching surface with boundary curve K corre- 
sponding to the boolean expression {ui V U2 V M3) A (ui V M2 V u^). 
The shaded surface F^-^ indicates the occurrence of ui in each of 
the clauses ci and C2. 

are identified with ki of the curves ci, . . . c„i on dF^, with one component of 
dFui identified with Cj for each occurrence of the literal Ui in the j*^ clause in C. 
Similarly, one curve of dF^^ is identified to the component of Fq labeled Mi, and 
the remaining ki boundary components are glued to ci, . . . Cm, with a component 
glued to Cj for each occurrence of the literal Ui in the j*'* clause of C. A total of 
three surface boundaries are identified along each of ui, . . . , u„, and exactly four 
surface boundaries are identified along each of ci, . . . , Cm, as in Figure 

Lemma 6. There is a truth assignment for U such that each clause in C has 
exactly one true literal if and only if there is a surface S with connected boundary 
and genus at most m + n and a continuous map f : S ^ B such that f\Qs is a 
homeomorphism onto K . 

Proof. Suppose there is a truth assignment for U such that each clause in C con- 
tains exactly one true literal. Form a surface S inside B by taking the union of 
Ff) and either Fu^ if Ui is true, or Fu^ if Ui is false. Then exactly two boundary 
components will be identified along each of the boundary components of Fq other 
than K itself, and K becomes the boundary of the resulting embedded surface S. 
There is a contribution of one to the genus from each of the literals, since F^. and 
Fg. each have genus one, and a contribution of one to the genus from each handle 
formed when a boundary component of a surface F„. or F^^ is glued to Fq along a 
curve Cj. The genus of S is therefore equal to m -I- n. 

Now suppose there is a surface Si of genus < n + m mapped continuously 
into B that has a single boundary component mapped homeomorphically to the 
boundary curve K. We will show in this case that there is a surface S4 with the 
same boundary, consisting of certain pieces of B identified along their boundaries, 
having genus precisely n + and containing, along each curve Cj, exactly one of 
the three pieces of surface joining Fq . 

The map /i of Si into B may be quite complicated, winding back and forth 
across B, but by standard transversality arguments we can homotop fi so that it 
is a union of homeomorphisms of subsurfaces of Si mapped homeomorphically to 
one of the pieces Fq, F„. , Fg. forming B. More precisely, we can perturb Si by 
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a small homotopy so that its intersection with the boundary components of Fq is 
transverse, and pulls back to a collection of simple closed curves on Si. If any of 
these curves bounds a disk in Si , then the disk is mapped into some subsurface X of 
B while the boundary curve is mapped to dX. Since X is not itself a disk, the disk 
can be homotoped into dX, and we can therefore homotop 5*1 in a neighborhood of 
this disk to remove a component of Si D f~^{dFo). After repeating finitely many 
times, each component of the complement of Si D /^^{OFq) in Si has non-positive 
Euler characteristic. The image of 5*1 in each piece Fq, f„. , F^. has an algebraic 
degree, which is either even or odd. This degree equals the number of pre-images 
in Si of a generic point in the piece. The degree of the map from 5*1 on Fq is odd, 
since K is the boundary of 5*1 and therefore 5*1 maps an odd number of times to 
points near K. The sum of the degrees along each of the pieces meeting a curve 
Ui or Cj is even, since 5*1 has no boundary along these curves. In particular, for 
each 1 < i < n, exactly one of -F„. , F^. has odd multiplicity in Si. Form a new 
surface 5*2 by taking the union of Fq and each of the pieces F^^ , F^. which have odd 
multiplicity in 5*1. We will show that x{S2) > xi^i)- The surface S2 is obtained 
from 5*1 by a series of operations that either discard a subsurface with non-positive 
Euler characteristic or replace a subsurface of Si that maps with odd degree to 
some piece F^ with a subsurface mapping homeomorphically to F^ . The Euler 
characteristic of a discarded subsurface of Si is smaller or equal to that of the 
subsurface that replaces it, so x(S'2) > xi^i) in either case. 

The collection of pieces among {Fu^, Fa^} that are in the image of S2 can be 
attached to Fq along common boundary curves among {wi, . . . w„}, forming a con- 
nected surface S3 of genus n. Finally, a connected surface S'4 with boundary K is 
obtained by identifying pairs of curves in dS^ that are mapped to the boundary 
components ci , . . . . Each c; then has either two or four curves in dS^ mapped to 
it, so there are either one or two identifications made along each c;. The choice of 
which pairs to identify, in case there are four curves mapped to c;, is not important. 
The Euler characteristic of S'4 is the same as that of S3 and S2 , and therefore greater 
or equal to that of Si. So the genus of S'4 is at most that of Si, genus(S'4) < n + m. 
Each identification of a, pair of surfa,CG boundaries along ci, . . . Cm contributes one 
to the genus of S4. There is at least one such identification along each of the m 
curves {cj}, though there may be two if each of the three surfaces meeting Fq along 
Cj has multiplicity one in S2. So identification of curves along ci,. . .Cm adds at 
least m to the genus of S4, and it follows that genus(S4) > n + m. Since we have 
seen that genus (S4) < n + m, equality must hold. Equality holds when exactly 
one of the three surface pieces meeting Fq along Cj has odd multiplicity for each 
1 < j < m. We then assign the value "TRUE" to a literal Ui if F^ is used in S4, 
and the value "FALSE" to Ui if F^. is used in S'4. This gives a truth assignment to 
U in which each clause in C has exactly one true literal. □ 

To show that 3-MANIFOLD KNOT GENUS is NP-hard we reduce in poly- 
nomial time an instance of the NP-hard problem ONE-IN-THREE SAT to an 
instance of 3-MANIFOLD KNOT GENUS. 

Proof of Theorem [B Given an instance of ONE-IN-THREE SAT form B as 
in Lemma El Then there is a truth assignment for U such that each clause in C 
has exactly one true literal if and only if K is the boundary of a surface of genus 
m + n mapped continuously into B. Form a 3-manifold N by thickening B so that 
it is embedded inside a triangulated 3-dimensional manifold with boundary. The 
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thickening process replaces each subsurface forming i? by a product of a surface 
with an interval, and then glues these surfaces together along portions of their 
boundaries, as indicated in Figure |21 The curve K remains on the boundary of TV, 
and there is a projection map p : N ^ B which fixes K. Form a closed manifold M 
by doubling N along its boundary, namely by taking two copies of N and identifying 
them along their boundaries by the identity map. Then M admits an involution t 
that fixes K, and has quotient N. 

To find a triangulation of N, we first describe an explicit triangulation of the 
branching surface B. An orientable surface with boundary has a triangulation with 
one vertex on each boundary component and no vertices in the interior. The number 
of triangles is 4g + 5c — 4, where g is the genus and c is the number of boundary 
components. We choose such triangulations for each of the 2n + 1 subsurfaces in 
the branching surface, and we match them together along boundary components to 
get a triangulation of B. So B has a triangulation in which the number of triangles 
is linearly bounded in n + m. We thicken the surfaces F^^ and F^^ , by taking their 
product with an interval. We form a cell structure of the thickened surface by 
dividing the product into prisms, products of a triangle and an interval. For the 
thickening of Fq, we start by doing the same and then go on to divide each interval 
into five subintervals. The top, middle and bottom subintervals are identified with 
the intervals from the three thickened surfaces meeting each thickened boundary 
component of Fq corresponding to a curve c^. Only the top and bottom intervals 
are identified with thickened surfaces from the other boundary components. We can 
therefore now glue these thickened surfaces together to get a 3-cell structure on N. 
Each cell in this structure is a prism. We form a closed 3-manifold M by doubling 
N along its boundary, gluing two copies of N together along their boundaries to 
obtain a 3-manifold with no boundary. 

Finally, we stellar subdivide the cell structure to get a triangulation. Each prism 
is divided into 14 tetrahedra, by dividing each rectangular face into four triangles by 
coning to a vertex in the center of each such face, and then coning the 14 triangles of 
the boundary of the prism to a vertex added to its center. The number of resulting 
simplices in AI is linearly bounded by n + m. 

We now check that if K bounds a surface of genus g < m + n with interior 
in M\K then it bounds a surface of the same genus in B. Suppose that F is an 
embedded surface in M with boundary K. If F does not already lie in N , then 
perturb it slightly so that the interior of F meets dN transversely in a finite number 
of simple closed curves and arcs. Using the involution r, refiect the portion of F 
not in N into iV, forming an immersed surface F' lying in N and with the same 
boundary as F. The interior of F' remains disjoint from K since K is fixed by 
the involution. The projection p{F') is a surface of genus g mapped into B with 
boundary K. So if K bounds an embedded surface of genus at most m + n in M 
then it bounds a surface of genus at most m + n mapped into B. The converse was 
shown in the proof of Lemma |S1 It then follows from Lemma |S1 that K bounds an 
embedded surface of genus at most m + n in M if and only if C is satisfiable with 
each clause containing exactly one true literal. 

The construction oi B,N and M described above each requires only a linear 
number of steps in the size of the instance of ONE-IN-THREE SAT with which 
we started, so that the reduction requires polynomial time. □ 
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4. Orbits of interval isometries 

In this section we develop a combinatorial procedure that will allow us to count 
the components of a normal surface. The procedure computes the number of orbits 
of a collection of k isometrics between subintervals of an interval [1, iV] C Z in time 
polynomial in klogN. By lining up the intersections of a normal surface with the 
edges of a triangulation, we obtain such subintervals. Arcs of the normal surface 
on the faces of the triangulation give rise to correspondences of these intersection 
points which are subinterval isometries. We can then apply the algorithm developed 
here to count the number of components of a normal surface. 

Assume that we have a set of integers {1,2,..., A''} and a collection of bijections, 
gi : [ai,bi] ^ [ci,di], 1 < i < k, either increasing or decreasing, that are called 
pairings. If a pairing identifies two intervals [a, b] and [c, d\ by sending a to c and b 
to d, we call it an orientation preserving pairing, and if it sends a to d and b to c, 
we call it orientation reversing. If a < ewe refer to [a, b] as the domain and [a, d\ as 
the range of the pairing. We work only with integers, and use the term "connected 
interval" to refer to the integers in a connected real interval. The width of an 
interval [a, b] with integer endpoints is 6 — a + 1, the number of integers it contains. 
The width of a pairing is the width of its domain or range, w = b — a + l = d — c+l. 
If the pairing preserves orientation, its translation distance t measures how far it 
moves points, so t = c — a = d — b. We can compose two pairings if the range of the 
first lies in the domain or range of the second. The collection of pairings generates 
a pseudogroup, under the operations of composition where defined, inverses, and 
restriction to subintervals. 

The interval [1, . . . , A] is divided into equivalence classes by the action of the 
pairings, which are called orbits. We are interested in the orbit structure of the 
collection of pairings, since with appropriate interpretation an orbit corresponds to 
a connected component of a normal surface. We introduce several simplification 
processes on the set of pairings in order to analyze the structure of the set of orbits. 

We introduce some terminology to describe the behavior of pairings. An interval 
is called static if it is in neither the domain nor the range of any pairing, so that its 
points are identified to no other points by pairings. Given a collection of pairings 
acting on the integers [1, . . . , A] , a pairing is said to be maximal if its range contains 
both A' and the range of any other pairing containing A^. More precisely, define a 
linear order on pairings using the lexicographical order (rfj, — Cj, — a^, —orientation), 
so that the maximal pairing has the highest upper endpoint, and among those 
with that endpoint the widest range, and among those with that range the biggest 
translation distance (if orientable). Finally, we say that a pairing g : [a,b] [c,d\ 
is periodic with period t if it is orientation preserving with translation distance t 
and a<c = a + t<b+l, so there is no gap between the domain and range. The 
combined interval [a, d] is then called a periodic interval of period t. 

The following lemma describes the orbits of a periodic pairing. 

Lemma 7. A periodic pairing has t = c — a orbits on [a, d]. 

Proof. Each point in [a, b] greater or equal to c lies in the range of g and can be 
mapped to a smaller point in [a,b] by a power of g~^. So each orbit on [a,d\ has 
a representative in [a,c — 1]. Since the congruence class modulo t of a point is 
preserved by g, each of the t integers in [a, c — 1] lies in a distinct orbit. These 
points uniquely represent the t orbits. □ 



THE COMPUTATIONAL COMPLEXITY OF KNOT GENUS AND SPANNING AREA 11 



We now show how to merge two pairings with sufficient overlap into a single 
pairing with the same orbits. 

Lemma 8. Let Ri be a periodic interval with pairing gi of period ti. Suppose that 
there exists an orientation preserving pairing g2 with translation distance t2 and 
an interval Ji C Ri such that Ji has width ti and g2{Ji) C Then the orbits 

of gi U 52 on Ri are the same as those of a single periodic action on Ri of period 
GCD(ii,i2). 

Proof. Let Ji be an interval in Ri of width ti that is in the domain of g2 and 
which is paired by 32 to an interval in Ri. Each point in Ri has a unique orbit 
representative in Ji under the action of gi on Ri. The interval g2{Ji) hes in Ri 
by assumption. For x G Ji let f{x) be the unique point in Ji obtained by carrying 
(72(2;) back to a point in Ji by a power of gi or g^^ ■ The effect of f{x) on Ji is a shift 
of t2 (mod <i). The orbits of / on Ji divide the points of Ji into congruence classes 
modulo GCD{ii,t2)- Neither gi nor g2 change the congruence class of a point mod 
(GCD(ti, t2)), so a subinterval of width GCD(ii,i2) in ^1 contains exactly one 
representative of each orbit of i?i under the action of 51 U 52 ■ The same orbits arise 
from a periodic action of period GCD(ti, t2)- □ 
The following is a special case of Lemma |H1 that applies when both gi and 172 are 
periodic pairings. 

Lemma 9. Let Ri and R2 he overlapping periodic intervals, associated with pair- 
ings 51, 52 having periods ti and t2. Suppose that width(Ri n -R2) > ii + ^2- Then 
the orbits of gi U 52 on Ri U R2 are the same as those of a single periodic pairing 
on Ri U i?2 of period GCD(ti,i2)- 

Proof. The leftmost interval Ji of width ti in i?i ni?2 is translated by g2 a distance 
of t2 to the right, to an interval J2 which lies in Ri fl i?2- Lemma |H1 then states 
that the action of gi U (72 on Ri is the same as those of a single periodic pairing 
on Ri of period GCD(ii,t2)- By symmetry the same result holds for the action of 
gi U 52 on i?2 . Thus the orbits of gi U g2 on i?i U R2 are the same as those of a 
single periodic pairing on Ri U R2 of period GCD(fi, i2)- D 
We state a consequence in a form which will be convenient for our applications. 

Lemma 10. Let gi,g2 be periodic pairings with periods ii,t2 and let Ji, J, J2 be 
intervals with gi{Ji) = J and g2{J) = J2- Suppose that J is contained in the union 
Ji U J2. Then the hypothesis of Lemma\^is satisfied and the orbits of gi U 32 cire 
the same as those of a single periodic pairing of period GCD{ti,t2), acting on the 
union of the periodic intervals of 51 , 52 • 

Proof. Let Ji ~ [ai,bi],J = [a, 6] and J2 — [a2,62]- We have 02 > ai and 
J C Ji U J2. Since Ji U J2 is a connected interval, we have that 02 < bi + 1. 
Note that ti =6 — 61 = a — ai and t2 ~ b2 ^ b — a2 ~ a. Then ti + t2 = 
{b - bi) + (02 - a) = (fe - a + 1) + (a2 - 61 - 1). The width of J is (6 - a + 1) and 
we have seen above that (02 — foi — 1) < 0. So the width of J is at least ti + 12 and 
Lemma implies the conclusion of the lemma. □ 
The orbit counting algorithm applies a series of modifications to a collection of 
pairings. We now describe these modifications. 
• Periodic merger 

The periodic merger operation replaces gi and 172 by a single periodic 
action on i?i U R2 of period GCD(ii,i2), as in Lemma El 
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Contraction 

The operation of contmction is performed on a static interval [r, s] . We 
eliminate this interval, replace [1, N] by [1, N — {s — r + 1)], and alter each 
gj by replacing any point a; in a domain or range which lies entirely to 
the right of s hy x — {s — r + 1). (This operation will lead to a decrease 
of s — r + 1 in the number of orbits, since the eliminated points are each 
unique representatives of an orbit.) 
Trimming 

The trimming operation simplifies an orientation reversing pairing whose 
domain and range overlap. Suppose that g : [a, b] [c, d] is a pairing with 
g{a) — d,g{b) = c and b > c. Define a new pairing g' : [a, (a + d)/2) 
{{a + d)/2,d] by restricting the domain and range of g, and say that g' is 
obtained from g by trimming. The domain and range of a trimmed pairing 
are disjoint. 
Truncation 

If an interval lies in the domain and range of exactly one pairing, then 
the interval can be "peeled off" without changing the orbit structure, in 
a way we now describe. The algorithm applies this operation to strip off 
points from the right of the interval [1, A'^]. 

When there is a pairing g : [a, b] — > [c, N] and a value A'"' with c < 
A''' + l < N, such that all points in the interval [A''' + l, iV] are in the range of 
no pairing other than g, then we can perform an operation called truncation 
of g. Truncation shortens the interval [l,iV] to the interval [1,A'''], and 
similarly shortens the domain and range of g. If g is orientation preserving, 
pairings other than g are unchanged, while g is eliminated entirely if c = 
N' + 1, or replaced by a shortened pairing g' : [a, 5 — (A^ — N')] [c, N'] 
if c < A'^'. We can perform this operation even if the interval [A^' + 1, Af] 
intersects both the range and the domain of g. 

If g is orientation reversing, truncation is applied only when g has disjoint 
domain and range (i.e. after trimming), and [A^' + 1, A^] is contained in the 
range. Suppose g : [a, 6] [c, A^] is a pairing with g{a) = N,g{b) = c 
and b < c, that [N' + 1, A^] is disjoint from the domains and ranges of all 
pairings other then g and that c < A'^' + 1. If A'^' + 1 = c then we eliminate 
g. Otherwise replace g by a shortened orientation reversing pairing, formed 
by restricting its domain to [a + N — N',b] and its range to [c, A''']. 
Transmission 

Transmission is an operation in which two pairings are composed. A 
pairing gi used to shift down the domain and range of a second pairing g2 
as much as possible. This operation will allow us to shift pairings leftwards 
from the right end of the interval [1, A^'] and subsequently apply truncation. 

If gi is orientation reversing and has overlapping domain and range, 
then as a first step in transmission we trim gi . Now consider a pairing gi , 
either orientation preserving or orientation reversing, and a second pairing 
92 whose range is contained in the range of gi . If the domain of 52 is not 
contained in the range of gi, form the composite map g2 = gi^ o g2, where 
r = 1 if 51 is orientation reversing and otherwise r > 1 is the largest integer 
such that gi^^^{[c2, ^^2]) is contained in the range of gi. The domain of 93 is 
the same as that of ^'2 in this case. If the domain of (?2 is also contained in the 
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range of gi then form the composite map g'2 = °32 °5i : gi''i[o-2, ^2]) 

{[^2, ^2]), where r is as above, s = 1 if 51 is orientation reversing, and 
otherwise s > 1 is the largest integer such that g^'*~''^([a2, 62]) is contained 
in the range of gi. The domain of is then that of 52 shifted left by g^^. 
The process of replacing 52 by g^^ ° 92 ° gf is called a transmission of 52 
by 91- 

We now construct a sequence of pseudogroups of pairings, terminating with the 
trivial pseudogroup acting on the empty set. At each stage, each orbit is associated 
to a unique orbit in the previous set of pairings, though the number of orbits may 
decrease. A counter is kept at each stage that records the total decrease in the 
number of orbits. Since the final pseudogroup has no orbits, the final value of this 
counter gives the initial number of orbits. 

Lemma 11. A contraction decreases the number of orbits by the width of the con- 
tracted interval. Altering a collection of pairings by any of the operations of periodic 
merger, trimming, truncation and transmission preserves the number of orbits. 

Proof. We refer to the collection of pairings gi, g2, . . . gk as G and to the new 
collection of pairings produced by one of the operations as G". 

The effect of a contraction is to shorten the interval [1,-/V] by removing points 
which are fixed by the entire collection of pairings. The number of orbits removed 
equals to the width of the contracted interval. 

A periodic merger joins two periodic intervals i?i and R2 and their pairings into 
one. Lemma |31 shows that the orbits in Ri U i?2 of gi U 52 are the same as those 
produced by a single periodic pairing g' . Suppose that x and y are points in [1, N] 
in the same orbit of the action of G = {51,32: ■■■9k}- Then h ■ x = y, where h 
is some finite word in the elements of G and their inverses. Wherever a gi or 52 
occurs in this word we can replace it by a power of g' , since g' translates by an 
amount that divides the translation distance of gi and 52- So x and y are in the 
same orbit of G' — {g' , 53, ... 5^}. Conversely suppose that that x and y are in the 
same orbit under G' . Then h' -x = y, where h' is some finite word in 5', 53, 54 ... 5^. 
Lemma 13 implies that wherever a g' occurs in h' we can replace it by some word in 
gi and 52, since the orbits of g' and gi U 52 coincide on the periodic interval of g' . 
So periodic mergers preserve orbits. 

Suppose next that gi : [ai,bi] [ci,di] is an orientation reversing pairing with 
gi{ai) = di,g{bi) = Ci < hi, and that : [0^,(0^ + di)/2) {{oi + di)/2),di] is 
obtained by trimming gi. If x and y are in the same orbit of G = {51,52, • • ■5fe}, 
then there is a sequence of points x = xi,X2, . . ■ ,Xr = y where each xj is the image 
of Xj-i under a pairing in G. We can replace an occurrence of gi or g~^ by 5' if the 
point Xj-i is smaller than (a^ + di)/2 and by g' ^ otherwise. So x and y remain in 
the same orbit under the action of G'. Conversely if x and y are in the same orbit 
under G' then replacing each occurrence of 5' by 5 and of 5' ^ by g~^ gives a word 
in G taking x to y. We conclude that trimming preserves orbits. 

Next consider the effect of a truncation. Suppose that all points in [7V' + 1, A^] are 
in the range of exactly one pairing 5 : [a, b] — > [c, N], and we truncate, shortening to 
an interval [1, N'] and either eliminating 5 or shortening it to 5' : [a, b— (A^ — A^')] 
[c. A''] in the orientation preserving case, or to 5' : [a + A^ — A^', b] [c, N'] in the 
orientation reversing case. Suppose that y — h ■ x where is a reduced product of 
pairings (a product that does not contain a subproduct of the form gj ■ gj^)^ We 
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can assume that y > x, as otherwise we can replace h by its inverse. Let h' be 
obtained from h by replacing all occurrences of g with g' . The successive images of 
X under the subwords of h can never enter and leave the interval [TV' + 1, A^], since 
they can only enter it under the action of some positive power of g, and only leave it 
under a negative power of g. These don't occur in succession in a reduced product. 
The image of a point z under gj is unchanged unless gj = g, and g{z) > N' + 1. 
So the image of x under h is the same as its image under h' unless y > N' . In that 
case h is of the form h ~ g\hi where hi ■ x < N' , r > 1, and y — g^hi ■ x > N'. So 
two points in [1, N'] are in the same orbit under G' if and only if they are in the 
same orbit under G. But each point in [N' + 1, A^] is in the same G orbit as a point 
in [1, N']. It follows that the number of orbits is unchanged by truncation. 

Finally suppose that a transmission of gj by gi replaces the pairing gj in G by 
g'j — g^^ ° 9] ° gt- li y — h ■ X where h is a word in G, then y — h' ■ x where 
h' is obtained by replacing every occurrence of gj with g^ o gj o g,^^. Similarly if 
y = h' ■ X where h' is a word in G", then y — h ■ x where h is obtained by replacing 
every occurrence of g'j with g~^ o gj o g^ . So transmissions preserve the collection 
of orbits and don't change their number. □ 

We now describe an algorithm which uses these operations to count the number 
of orbits of a pseudogroup of pairings acting on [1, A'^]. We initially set a counter 
for the number of orbits to zero. A pairing gi is said to be maximal if di — N 
and [ci,N] contains the range of any other pairing with an endpoint at N. Let 
S = {1, 2, . . . N} and let gi : [a^, bi] — > [ci, di], 1 < i < k he a collection of pairings 
between subintervals of S. The algorithm will contract S and reduce the number 
of pairings. It keeps a running count of the number of orbits detected in an integer 
referred to as the orbit counter. Denote by N' the current size of the interval as we 
proceed. The algorithm repeats the following steps, reducing N' and k, until there 
are no points remaining. 

Orbit counting algorithm: 

(1) Delete any pairings that are restrictions of the identity. 

(2) Make any possible contractions and, if any exist, increment the orbit counter 
by the sum of the number of points deleted by the contractions. If the num- 
ber of pairings remaining is zero, output the number of orbits and stop. 

(3) Trim all orientation reversing pairings whose domain and range overlap. 

(4) Search for pairs of periodic pairings gi and gj whose domains and ranges 
satisfy the condition of Lemma El If any such pair exists then perform a 
merger as in Lemma El replacing gi and gj by a single periodic pairing, 
with translation distance GCD (ti,tj). The new pairing acts on the union 
of the domains and ranges of gi and gj. Repeat until no mergers can be 
performed. 

(5) Find a maximal gi. For each gj ^ gi whose range is contained in [ci,N'], 
transmit gj by gi. 

(6) Find the smallest value of c such that the interval [c, N'] intersects the 
range of exactly one pairing. Truncate the pairing whose range contains 
the interval [c,N']. 

Theorem 12. The orbit counting algorithm gives the number of orbits of the action 
of the pairings {gi}i^i on [1, . . . , N] in time bounded by a polynomial in kiogN. 
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Proof. We first check that the orbit counting algorithm correctly counts the num- 
ber of orbits. In Step (1), deleting a pairing which is the identity on its domain does 
not change the number of orbits. In Step (2), contracting a static interval removes 
a number of points that are unique orbit representatives, and the count of these is 
added to the running total kept in the orbit counter. If all the points have been 
removed, then there are no more orbits to count and the orbit count is complete. 
Lemma 1111 shows that the operations of transmission, trimming, truncation and 
merger occurring in Steps (3), (5) and (6) do not change the number of orbits of the 
collection of pairings. 

Mergers carried out in Step (4) occur when the conditions of Lemma |31 are 
satisfied, and Lemma^lshows that they preserve the orbit structure of the collection 
of pairings. A merger reduces the number of pairings by one, replacing two pairings 
gi and gj by a single periodic pairing acting on the union of their periodic intervals. 

In each cycle through these steps the interval width decreases by at least one, 
either in Step (2) or in Step(6). It follows that the algorithm terminates, yielding 
a count of the number of orbits after a number of steps bounded by the width N of 
the interval. We will obtain a much better bound. To do so, we define a complexity 
which decreases as we iterate the above steps. Recall that Wi = bi — ai + l, 1 < i < k 
and k is the number of pairings. The complexity X is defined to be 

A; 

1=1 

The process of executing Steps (l)-(7) in turn is called a cycle. We will show that 
when we run through 2k cycles, X is reduced by a factor of at least two. See the 
remark at the end of this section for a geometric interpretation of this complexity. 

Call an interval [x,y] a Z -close interval if it is the domain or range of a pairing 
of subintervals of [1, Z] and ify — a; + l > {Z — x + l)/2, or equivalently, if y > 
Z/2 + x/2 — l/2. Being Z-close corresponds to being relatively close to Z. It means 
that there is no room for another interval of the same size to the right of y. The 
value of Z is initially set to N, and as the algorithm proceeds it is reset to the 
current interval width N' each time the number of pairings k is decreased. The 
number of pairings decreases when a merger occurs, or when a pairing is truncated 
to zero width, or when a pairing is transmitted to become trivial (a restriction of 
the identity) and then eliminated. 

Claim 1. The union of two Z-close intervals of equal width is a connected interval. 

Proof. Suppose that [a, b] and [c, d] are equal width Z-close intervals with c > a. 
Then b-a = d-c, b-a + l> (Z - a + l)/2 and d - c + I > {Z - c + l)/2, so 
b > Z/2 + a/2 - 1/2 > d/2 + a/2 - 1/2 = 6/2 + c/2 - 1/2 implying that 6 > c - 1 
and [a, b] U [c, d] is connected. □ 

Claim 2. Suppose that the domain [a, b] C [1, Z] of a pairing g is not Z-close. 
Then the image [x',y'] of an interval [x^y] under g~^ is not Z-close. 

Proof. Suppose that [a;', y'] is the image of [x, y] under g~^ , where g : [a, b] — > [c, d]. 
By assumption we have that b < Z/2 -\- a/2 — 1/2. The interval [x' , y'] is contained 
in [a, b] so a < x' < y' < b. Then y' < Z/2 + a/2 - 1/2 < Z/2 + x' /2 - 1/2 and 
[x',y'] is not Z-close. □ 
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Claim 3. After a series of five cycles either the number of Z -close intervals de- 
creases or the number of pairings decreases. 

Proof. In each cycle at least one inaximal pairing is truncated. Suppose that 
during a series of five cycles, five successive maximal pairings gi, 92, 93, 94, 95 occur 
in turn in the orbit counting algorithm. The initial maximal pairing gi is truncated 
in the first cycle. Eventually gi stops being maximal, and it is then transmitted by 
the new maximal pairing g2 in the next cycle. The cycle 52 is then itself truncated 
until no longer maximal and it in turn is then transmitted by 173 and so on. When 
(7i+i transmits gi, the pairing of [0^,6^] with [a, di] is replaced by a pairing of [a[, b'^ 
with [c^,ci-] and one of the following three cases occurs. 

(1) The range [q, di] of gi is transmitted to a non Z-close interval by gi+i- 

(2) The range [ci,di] is transmitted to a .^-close interval [cj,rfi] by gi+i, and 
the domain [a^, b'j] of the transmitted pairing is Z-close. 

(3) The range [ci,di] is transmitted to a Z-close interval [c^,(i^] by gi+i, and 
the domain [a'^, b'^\ of the transmitted pairing is not .Z-close. 

We first consider the case where each of gi, g2,andg^ is orientation preserving. 
Setting Z to the initial interval width N' and noting that gi is initially maximal, we 
see that the range of gi is initially Z-close, with Z = N' = di. Truncation reduces 
the range to of gi to [01,^2], at which point the pairing 32 becomes maximal. We 
can assume the range [ci, ^2] of gi is still Z-close, or we are done. In the next cycle 
the interval [ci,c?2] is transmitted by 32 to an interval [c']^,d']^]. If Case (1) applies, 
this new interval is not Z-close and the number of Z-close intervals has decreased. If 
Case (2), the three intervals [ci, ^2], [c^ , d'^] and [a'^, b'^] are all Z-close and of equal 
width, and therefore satisfy the hypothesis of Claim Lemma El then implies 
that two pairings can be merged, and the number of pairings decreases during the 
execution of Step (4) in the next cycle. If Case (3) occurs, the maximal pairing 172 
is truncated in the second cycle. We can assume its domain remains Z-close, or 
we are done. When the next pairing g^ becomes maximal and transmits 52 in the 
third cycle, we must fall into one of the first two cases. So the number of Z-close 
intervals is decreased in the third cycle or the number of pairings decreases in a 
merger during a subsequent application of Step (4) in the fourth cycle. 

Now we consider the possibility that one of 52, andg^ is orientation reversing. 
The domain and range of a trimmed orientation reversing pairing are disjoint, and 
hence any interval transmitted by an orientation reversing pairing is not Z-close. If 
any of g2, gsorg^ are orientation reversing then it transmits the previously maximal 
pairing's range to an interval that is not Z-close and the number of Z-close pairings 
decreases. If non is orientation reversing then there is a sequence of three successive 
orientation preserving pairings and the previous argument applies. □ 

Claim 4. Suppose that Z is set to the current interval size N' and that after a 
series of truncations in which the maximal pairings are successively 31,92, ■ ■ • ,gr, 
no Z-close intervals remain. Then the complexity X is reduced by a factor of at 
least two. 

Proof. Truncation of a maximal pairing gi results in its range [ci,di], of width Wi, 
being reduced to a shorter interval [ci, c?'J of width w[. This reduces the complexity 
X by a, factor of 

■14 _ d[ -c, + l 
Wi di - Ci + \' 
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The maximal pairing changes after a truncation if the range of g[ has truncated 
sufficiently so that it is contained in the range of gi+i, so that Ci+i < c[ and 
di+i = d\. Define d" by 

= Ci - f + [w,+i)- — . 

Z - Ci+i + f 

and define 

w'l = d'l - c, + 1. 
Differentiation shows that the function 

J., . X - Cj + 1 



X - Ci+l + f 

is increasing with x for x> Ci, since q > Ci+i. So 

Wi = dj - c,j + f = dj+i - Cj + f = (wj+i) 

Wi+l 

. , [di+i - Cj + f ) Z - Ci + 1 

= {Wi+i)— — < (Wi+l)- 



(di+i - Cj+i + f ) Z - Ci+i + 1 

= d'l -c, + l^ w'l 

After a series of truncations of the pairings 51, 172, • ■ ■ , 5r, the complexity X is 
multiplied by a factor of 

WlJ \W2j \Wr 



< 



w'4 \ I w'' 



Wl J \W2 J \Wr-l J \Wr 

d" - Ci + 1 \ / 4' - C2 + 1 \ / 4-1 - Cr_i + 1 \ / w'j. 



Wl J \ W2 J \ Wr-l / V'"^ 

W2\ f Z — Ci + \ \ ( W,. \ f Z — Cr-1 + 1 \ f W. 



Wl/ — C2 + I/ \Wr-l J \ Z — Cr + 1 J \Wr 

' Z-Ci + 1 \ 
Z - Cr + I J \Wi 

Now Z = di, since gi was the first maximal pairing truncated, so Z — ci + 1 = 
di — ci -\- 1 — Wl. By assumption, there are no Z-close intervals remaining after 
gr is truncated, so that [cr,c?5,] is not Z-close. We then have by the definition of 
Z-close that 

^-ci + l \ fK\ ^ < ^ d', - c, + 1 ^ ^ 

Z — Cr + 1 J \Wi J Z — Cr + I Z ~ Cr + 1 

It follows that X is multiplied by a factor smaller than 1/2. □ 
Proof of Theorem 1121 The operations involved in one cycle consist of compar- 
isons, additions, subtractions and computing greatest common divisors of the k 
pairings of G. The number of these operations occurring in each cycle is linear 
in fc^logiV. Pairings are described by pairs of intervals, whose boundary points 
are integers of size at most N. So the running time of each cycle is polynomial in 
A: log TV. 

If the number of pairings k decreases during a cycle, then a pairing has been 
eliminated because its width has truncated to zero, because it has been transmitted 
and become the identity on its domain, or because two pairings have merged. When 
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two pairings 51,(72 merge, the value of k decreases by one. The product (wi)(w2) 
which occurs in X is replaced by the width of the new pairing, which is at most 
{wi + W2) < 2wiW2- Therefore X — 4'^ J|*L]^ Wi decreases by a factor of at least 
two after a merger. 

Since each pairing has width at most N, the initial complexity is bounded above 
by N^A''. By Claim O each time we run through five cycles either the number of 
pairings decreases or the number of Z-close intervals decreases. There are at most 
k Z-close intervals initially, so after 5k cycles either there is a reduction in the 
number of pairings or there are no Z-close intervals remaining. In the second case 
the complexity X has decreased by a factor of at least two, by Claim 01 

So X decreases by a factor of at least two after 5k cycles, and the complexity 
reduces to zero after at most fc(2 + log2 N) successive series of 5k cycles, or after 
at most 5fc^(2 + logj iV) cycles. Each cycle runs in time polynomial in fclogA^, so 
the total running time is also polynomial in fclogiV. □ 

We now apply Theorem ll2l to count the number of components of a normal curve 
or normal surface. An obvious algorithm to count components proceeds by marking 
vertices connected by common edges until all vertices in a component are reached. 
This procedure takes time linear in the number of edges of the curve. This is equal 
to the sum of the normal coordinates (called W below), but we can achieve an 
exponential improvement. We first look at normal curves. 

Corollary 13. Let F be a surface with a triangulation T containing t triangles and 
let J be a normal curve in F with normal coordinates summing to W . There is a 
procedure for counting the number of components of 7 that runs in time polynomial 
in t log W . 

Proof. The 1-skeleton of T contains e edges, where e < 3t. Fix once and for all 
an ordering of these edges. A normal curve 7 intersects each edge of T in a finite 
number of points. Set N to be the weight W of 7, the sum of the number of 
intersection points of 7 with all the edges of T. Label the intersections of 7 and the 
first edge of the 1-skeleton by the integers 1, 2, . . . , ii, the intersections of 7 and the 
jth edge hy ij-i + l,ij-i + 2, ... ,ij, and the intersection of 7 and the e*'' edge of T 
by «e-i + 1, ie~i + 2, . . . , Zg. Then ig = N. Each triangular face of T has three sets 
of arcs pairing points of [1, A'^], with one set running between each pair of edges of 
the face. To each set of arcs we associate a pairing between the intervals at either 
end of the arcs, as in Figure 0] All the pairings are orientation reversing in this 
example. In general some will be orientation preserving, as the edge orientations 
on any triangle can be arbitrary. 

The number of connected components of 7 is the same as the number of orbits 
of an action of a collection of pairings on 7 H T^-^\ where two points are paired if 
they are connected by an edge of 7 lying in a triangle. This is precisely the number 
returned by the orbit counting algorithm. The number of pairings is at most 3i 
and the interval size is W. Applying Theorem ll2l we can determine the number of 
components of the normal curve 7 in time polynomial in tlogW. □ 

A similar argument applies to normal surfaces. 

Corollary 14. Let M be a 3-manifold with a triangulation T containing t tetra- 
hedra and let F be a normal surface in M whose normal coordinates sum to W . 
There is a procedure for counting the number of components of F which runs in 
time polynomial in tlogW . 
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Figure 4. Normal arcs on a triangular face give three pairings of 

ordered intervals: gi : [rii + ri3, 712] [^3,713 + ri2 - 1], 52 : [ns + 
r 12, tia] [n5,n5+r23-l], and^a : [n5+r23,n6] [ni, ni+ri3-l] 
The normal curve coordinates associated to this face are ri2,ri3 
and r23. In a normal surface, ri2,ri3 and r23 are each a sum of 
two normal coordinates. 

Proof. The 2-skeleton of T contains at most At faces and the 1-skeleton contains 

e < 6t edges. Set to equal the weight W of F, the total number of points in 
which it intersects the 1-skeleton. Order the edges of T^^^ in an arbitrary way and 
label the intersections of F and the jth edge of T^^^ by + 1, ij-i + 2,...,ij. 
Again we have ie = N. To a pair of edges on a triangular face of the 2-skclcton 
we associate a pairing between the intervals at either end of the corresponding set 
of arcs. There are at most three pairings for each face, and the number of faces is 
at most 4t, so the number of pairings is bounded above by 12t. These pairings are 
determined by the normal coordinates of F. 

For a normal surface, the number of connected components of F is the same 
as the number of components of _F n T^^'' , since every component of F intersects 
the 1-skeleton of T, and if two points on can be joined by a path 

then that path can be homotoped into the 2-skeleton of the triangulation. So the 
the number of components equals the number of orbits of _F n T*-^^ under pairings 
that identify two points connected by an edge of F n T^^' contained in a face of 
the 2-skeleton. This number is precisely what is computed by the orbit counting 
algorithm. Therefore we can determine the number of components of the normal 
surface F in time polynomial in tlogW. □ 
Remark: There is a motivating geometrical construction behind the combinatorics 
of the orbit counting algorithm which we informally explain. We can associate to 
each pairing gi : [ai,6i] —^ [ci,di] a pair of "transmission towers" in the upper half- 
plane, with one tower being the vertical line segment from (ai,0) to {ai,Wi) and 
the other the vertical segment from (cj,0) to {ci,Wi). These towers capture the 
information contained in the pairings. Points in the domain and range of a pairing 
beam up to the tower leftward at a 45 degree angle, then beam to the paired tower 
(either straight across, if the pairing preserves orientation, or crossing if not) then 
down again. Assign a "cost" to each transmission tower equal to the hyperbolic 
length in the upper half space model from y = 1/e to the top of the tower. This 
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cost equals log(di — a + 1) + 1, and the sum of all these costs is essentially the 
logarithm of the complexity X used above. 

The counting algorithm starts with a Euclidean line emerging from the right 
endpoint of a Z-close interval and going upward at a 45 degree angle. In hyperbolic 
terms, this is an equidistant curve from the geodesic from the right endpoint to the 
point at infinity. Initially, we may assume that this sweep line hits the top of at 
least one transmission tower. We make the highest of these the Z-close tower, and 
use it to beam all the other towers to its paired interval. If the two domains overlap, 
we use the highest power of the transmission that can be applied. Eventually the 
triangle between the current tower and the sweep line is vacant. We sweep leftward 
with this equidistant curve. As long as the sweep line hits only one tower, we 
truncate it. 

Consider a second equidistant curve P through a point x on R, with slope -1/2. 
The hyperbolic distance between the slope -1 and slope -1/2 lines is a constant 
equal to about 0.49. Call the region between these two lines a zone. If two paired 
towers have their tops in a single zone, their domains overlap. If no towers are 
completely truncated away between the time when the sweep line hits a point x 
and when there are no towers intersecting the zone above P, each equidistant curve 
between the two given ones hits a transmission tower based at a point that will be 
removed by truncation. So in this time either the number of towers has decreased 
or the sum of all the tower costs has decreased by at least 0.49. 

If we merge transmission towers when possible, a complexity based on the cost 
decreases sufficiently fast to give a polynomial time algorithm. The calculations in 
Theorem 1121 implemented this geometric picture. 

5. 3-MANIFOLD KNOT GENUS is NP. 

In this section we establish that 3-MANIFOLD KNOT GENUS is NP. 
Proof of Theorem |2j We begin with a simplicial complex consisting of t tetrahe- 
dra whose faces are identified in pairs, and a collection K of edges in the 1-skeleton 
of this complex. While there are alternate formats in which a knot and a 3-manifold 
may be presented, all reasonable ones appear to be transformable to one another 
in polynomial time. 

In time polynomial in t we can check that the link of each vertex is connected 
and has Euler characteristic two, which means that it is a sphere, and that the link 
of each edge is a connected curve. These are the necessary conditions to ensure 
that the underlying space of the complex is a 3-manifold M. Similarly we can 
check in time polynomial in t that the edges of K form a simple closed curve in 
M, and that this curve represents a trivial element of the first homology group of 
M with integer coefficients. We then form the second barycentric subdivision of 
the triangulation of M, replacing each tetrahedron by 576 tetrahedra. Removing 
all closed tetrahedra that meet K results in a 3-manifold Mk with a single torus 
boundary component, the "peripheral torus" that surrounds the knot K. The knot 
K bounds a surface of genus 5 in M if and only if there is a surface in Mk of 
genus g with a single boundary component that is an essential curve on BMk- We 
restate our problem in the triangulated manifold Mk as the question of whether 
there exists in Mk an orientable surface of genus g with a single essential boundary 
component on dMx- By Schubert if such a surface exists, then such a surface 
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exists among the fundamental normal surfaces in Mk- The certificate consists of 
an integer vector w in 17^ giving the normal surface coordinates of this surface. 

Recall that not all vectors in correspond to normal surfaces. It is necessary 
that w satisfies the matching, positivity and quadrilateral conditions. These con- 
ditions can be checked in time which is linear in t. So we can verify that there is a 
normal surface F with w — v(i^). 

To verify that dF is essential on the peripheral torus, we include in the certificate 
a non-trivial cycle in the 1-skeleton of OMk that intersects dF in an odd number of 
points. Such a cycle can be found in the 1-skelcton of dMx, since curves embedded 
in the 1-skeleton generate its first homology. Odd intersections with such a cycle 
implies that dF is non-separating on dMx, and in particular does not bound a disk 
on OMk, ensuring that dF corresponds to a longitude curve parallel to K in Nk- 
Using the orbit counting algorithm, we can count the number of components of 
the normal surface F and verify that F is connected in time polynomial in t log W, 
where W is the weight of F the number of points in which F meets the 1-skeleton 
of Mk- We apply Corollary 1141 to verify that dF is connected in time bounded by 
a polynomial in tlogW. (This last step can be avoided. The number of essential 
boundary components must be odd, since dF represents a non-trivial element in 
Hi^dMx] Z2), and an even number of them can be removed by joining adjacent 
pairs of curves with annuli on dMx- This gives a surface with one boundary 
component and the same Euler characteristic. Inessential boundary curves can be 
capped while increasing Euler characteristic, which gives lower genus.) 

To check that F is orientable, we take the vector 2v(F) that doubles each coordi- 
nate of the normal surface F , and apply the orbit counting algorithm to determine 
if the corresponding normal surface F is connected. Since M is orientable, F is 
connected if F is connected and non-orientable, and has two components if F is 
connected and orientable. Thus we can verify if F is orientable in time which is 
linear in t log W . 

The Euler Characteristic x(-FO is determined by the number of vertices, edges 
and faces of _F, which are computable from its normal coordinates in time which is 
linear in t. Following |11| . we let ti be the number of tetrahedra containing edge 
and set e.y — 1 if the edge meets the j*'* normal disk. Then the normal surface 
F with coordinates {vj} has x{F) = i.^/'^)h - <^iF) + wt{F) where a{F) = J^^j 
and wt{F) — J2i j ^ij^i/U. The values of eij,ti are determined by the triangulation 
and are independent of F. They can be computed in time polynomial in t. Since 
_F is a connected orientable surface with one boundary component the genus of F 
is (1 — x)/2- Thus we can determine the genus of F in time which is linear in t. 

Theorem 0] implies that the normal coordinates of this surface are at most i2^*+^. 
There are 7t normal coordinates, and each represents a triangle or quadrilateral, so 
that the total number of intersections with the 1-skeleton satisfies N < 28i^2''*~'"^. 
In particular, log is bounded above by a polynomial in t. So the fact that _F is a 
spanning surface for K can be verified in time polynomial in t. □ 

6. An EXTENDED COUNTING ALGORITHM 

In this section we develop a generalized version of the orbit counting algorithm, 
that counts not only the number of orbits of a collection of isometrics between 
subintervals of an interval, but also more general quantities which are useful in ap- 
plications. For example, we can use the extended algorithm to answer the following 
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question: Given a normal surface and a triangulation, how many times does each 
component of the surface intersect a fixed edged of the triangulation? The extended 
algorithm allows one to effectively compute the normal coordinates and the Euler 
characteristic of each connected component of the surface, hence the genus, even 
when there are exponentially many components. To carry out such computations, 
we extend the previous analysis to pairings of weighted intervals, in which each 
point of the interval has associated to it a vector in 1'^. We are interested in the 
sum of these vectors over an orbit, the orbit weights. 

Consider again a pseudogroup of interval isometrics acting on [1, . . . , iV] . We will 
assume that there is given as input a nonnegative weight function 2; : [1, . . . , TV] ^ 
Z+d, associating to each element of [1, . . . , iV] a vector in Z^J. satisfying the following 
condition: the weight at successive points j and j + 1 changes at most 4/c times. 
In our application, k will be the initial number of pairings, and 4fc gives an upper 
bound on how many times an endpoint of a domain or range of a pairing is reached 
as one moves across \\.,. . . ,N\. The algorithm proceeds as before, while maintaining 
data on the orbit weights. We keep track of the orbit weights by maintaining two 
lists of weighted subintervals. The first 

-^ = {([Pl,9l],2:i),---,(bm,gm],^;m)}, 1 < Pi < % < 

records the current weight values at each point in [1, AT], with initially equal to 
N. This list is updated as the algorithm proceeds. Points in the interval 
have constant weight Zj £ Z'^, and there are at most Ak such intervals. A second 
list of t subintervals 

L' = {{[n, si],vi), . . . , {[n, st\,vt)}, 1 < n < Si < n+i 

consists of a collection of intervals [r^ , Si] paired with a vector Vi S . This pair 
represents (sj — + 1) orbits, one for each point in the interval [rj, Sj], and to each 
of these orbits is assigned the orbit weight Uj. Initially empty, at the algorithm's 
conclusion L' records the total number of orbits ^^(si — + 1), along with the 

i 

orbit weight assigned to each of the (sj — rj + 1) orbits in the interval [rj, Sj]. 

Wc define an additional operation, called transferring weights by a pairing g. 
Suppose that g : [a, b] — > [c, d] is a pairing and that [c, d] carries n different weights, 
given by the list 

{([C = ri,Si],Vi),...,{[rn,Sn = d],Vn)}, l<ri<Si< n+i. 

The weight function can be split into constant functions on n subintervals of [c, d] , 

where l<n<d — c+1. 

The transfer operation sets the weights on [c,d] to zero and keeps the orbit 
weights the same by translating the weight vectors of [c, rf] to smaller orbit repre- 
sentatives, as below: 

Case 1: c/ is orientation preserving and b < c. Set the weights on [c, d] to zero 
and for each 1 < j < n, add vj to g~^{[rj,Sj]). 

Case 2: g is orientation preserving and b > c. Then g : [a, b] [c, d] is a periodic 
pairing of period t = c — a. Wc set the weights of points in [c, d] to zero and adjust 
weights in [a, c — 1] to preserve the orbit weight. The points in [c, d\ have weights 
given by the intersection of [c, d] with intervals in L. These weights are described 
by 

{([c = ri,si],vi), ([r„,s„ = d],Vn)}, 1 < n < Si < n+i. 



THE COMPUTATIONAL COMPLEXITY OF KNOT GENUS AND SPANNING AREA 23 



For each interval {[rj, Sj],Vj), 1 < j < n, of width Wj and constant weight Vj, add 
[[— to the weight of each point in [a, c — 1], and add an additional Vj to the 

weight of each point in [a, c— 1] that is congruent mod(t) to a point in Sj], 
if any such point exists. 

Case 3: g is orientation reversing. Wo first trim g. This docs not affect the orbits 
or the orbit weights. For a trimmed, orientation reversing pairing g : [a,b] [c, d], 
set the weights on [c, d\ to zero and for each 1 < j < n add Vj to g~^{[rj,Sj]). 

Lemma 15. The operation of transfer sets the weights on [c, d] to zero and pre- 
serves the orbit weights of a collection of pairings. The number of distinct weights 
taken by the weight function on [1,N] increases by at most four following a transfer 
operation. 

Proof. In Cases (1) and (3), the decrease in the weight function at one point x in 

an orbit is exactly offset by an equal increase at the point g~^{x) in the same orbit. 

In Case (2), each point in [a,d] is in the orbit of a unique point in [a,t — 1] = 
[a, c— 1] under the iterates of g. The total weight of an orbit within a periodic pairing 
is transferred to the orbit representative in this initial subinterval by applying 
powers of g. Adding weight vj to an orbit representative in [a, c — 1] of each point 
in [rj,Sj while setting the weight at that point to zero, preserves the orbit weight. 
The resulting weights on [a, c — 1] are gotten by adding appropriate multiples of Vj, 
with the factor being the number of orbits of a point in [a, c — 1] that lie in [rj, Sj]. 

The number of orbits under of a point in [a, c — 1] that lie in [r^, Sj] is [[~^]] '^^ 

one more than this for points whose orbit hits the last wj {mod t) points of [rj, Sj]. 
It follows that the transfer operation in Case (2) preserves the number of orbits 
and the orbit weight. 

Setting the weights on [c, d] to zero can cause at most two new points where 
a weight change occurs. The transferred weights from a constant weight interval 
[rj,Sj] result in a net increase of at most two pairs of successive points where 
the weight changes, at the preimages of its two endpoints. In Cases (1) and (3), 
transferred weights from an interval with non-constant weights rcsiilts in an increase 
in the number of weight changes in the domain of g. The increase in the number 
of weight changes in the domain is exactly canceled by the decrease in the number 
of weight changes in the range of g, except possibly for two extra weight changes 
at the boundary points a, b of the domain. In Case (2) the same holds, but with 
[a, c — 1] replacing the domain. In each case the number of constant weight intervals 
m' is increased by at most four during a transfer operation. □ 

We now describe the modified algorithm. Again N' represents the current inter- 
val length, and we set m' to be the current number of constant weight intervals. 

Weighted orbit counting algorithm: Let {gi : [aj,6j] — > [cj,dj], 1 < i < k} he 
a collection of pairings between subintervals of {1, 2, . . . N}, and let 

L = {{\pi,qi], Zi), . . . ,{\pm,qm\, Zm)}, ^ < Pi < <Pi+l <N 

be a list representing a collection of weights on [1, A/'], with the weight on \pj,qj] 
equal to Zj G Z''. Initialize a second weight list V to be empty. The algorithm 
proceeds as before, reducing the interval size A^ until it reaches zero, but this time 
keeping track of orbit weights by maintaining the lists L,L' . 
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(1) Search through the pairings and delete any pairings which are restrictions 
of the identity. Leave the weight hsts L, L' unchanged. 

(2) Search for and contract static intervals. If the interval [r, s] is contracted, 
has constant weight z, and is not contained in a larger contracted interval 
with the same weight z, add an interval {[N' + 1, N' + s — r +1], z) of width 
s — r + 1 to the end of L' , with associated weight z. Alter L by replacing 
[1, N'] by [1, N' — {s ~ r + 1)], and altering each gi by replacing any point 
X in the domain or range of gi with x > s by x — (s — r + 1). Replace the 
weight function by a new weight function w' , which at points x > s satisfies 
w'{x) — w{x + {s — r + 1)) and agrees with w at points x < r. 

(3) Trim all orientation reversing pairings whose domain and range overlap. 
Leave the weight lists L, L' unchanged. 

(4) Search for pairs of periodic pairings gi and gj whose domains and ranges 
satisfy the condition of Lemma |3 If any such pair exists then perform a 
merger as in Lemma |51 replacing gi and gj by a single periodic pairing, 
with translation distance GCD(ii,ij), acting on the union of the domains 
and ranges of gi and gj. Leave the weight lists L,L' unchanged. Repeat 
until not mergers can be performed. 

(5) Find a maximal gi. For each gj with j ^ i, if the range of gj is contained 
in [ci,N'], transmit gj by gi. Leave the weight lists L,L' unchanged. 

(6) Find the smallest value of c such that the interval [c, N'] intersects the 
range of at most one pairing gi, with gi : [ai,bi] — > [ci,N']. Transfer the 
weights on [q, N'] by gi, and then truncate the pairing gi. 

(7) If the interval size N' has decreased to zero, output the list L' and stop. 
Otherwise start again with Step (1). 

For a Z''- valued function on [1, 2, . . . , A^] whose values are given by the list L = 
{{[pi,qi],zi), . . . , {[p„i,q„i],z„r)}, I < Pi < Qi < Pi+i, define the total weight of L 
to be J2Zi l^'l- 

Theorem 16. Suppose there is a pseudogroup generated by k pairings with Un- 
valued weights, {gi}i^i on [1,N] such that there is a partition of [1,N] into m 
disjoint subintervals in which the weights are constant, and such that the total 
weight is at most D. Then the weighted orbit counting algorithm outputs a list 
with one point for each orbit and corresponding orbit weights, and runs in time 
polynomial in kmd\ogD\ogN. 

Proof. We will check that the running time of the algorithm is larger than that of 
the previous unweighted version by a factor which is a polynomial in mdlogD. 

The proof that the algorithm terminates is the same as that given for Theorem ll2l 
There is some extra overhead involved in keeping track of weights that modifies the 
calculation of the running time. We indicate these additional calculations below. 
We now check that at each step in the algorithm the orbit weight is unchanged 
for any orbit remaining in L, and that eliminated orbits have their orbit weights 
correctly recorded in L'. 

As we run through the steps of the algorithm. Steps (1),(3),(4),(5) and (7) pre- 
serve the orbit structure, the weight function and the interval [1,A^'], so neither 
of the lists L, L' are changed. The number of constant weight sub-intervals is also 
unchanged. 

Step (2), contraction, does change the orbit structure, and also shortens [1, A^']. 
In Step (2) the procedure adds the eliminated orbits and their weights to L' . The 
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number of constant weight intervals is not increased, and may be decreased. Main- 
taining the two hsts requires at most 0{md\ogD) additional steps. 

In Step (6), truncation, points are eliminated from the end of the interval [1, N']. 
However since we first transfer the weights of these points, the eliminated points 
all have weight zero and the weight of an orbit is unaffected. 

The number of steps involved in resetting the weights in L for a transfer operation 
is given by a polynomial in m'dlogDlogN. Since m' increases by at most four 
at each of the polynomially many steps of the algorithm, m' is bounded by a 
polynomial in mklogN. 

Combining the running time of each of the steps, whose number is given by a 
polynomial in fclogA'^, gives a polynomial in kmdlogDlogN for the total running 
time. □ 

Corollary 17. Let M be a 3-manifold with a triangulation T containing t tetra- 
hedra and let F be a normal surface in M of total weight W . There is a procedure 
for counting the number of components of F and determining the topology of each 
component which runs in time polynomial in tlogW. 

Proof. We begin as in Corollarvll4lbv assigning an integer in [1, N] to each point 
of intersection between the normal surface and an edge of the triangulation, where 
N is the total number of intersections of F with the 1-skeleton, and again associate 
three pairings to each face of the triangulation, one to each pair of edges in the 
face. The number of pairings that results is bounded above by I2t. 

We next define a weight function w{x) which assigns integer weights (zi, Z2, . . . , 274) 
to each point in [1, A^]. Initially Zi is set to zero for all i at all points a; G [l, N]. A 
tetrahedron can have as many as five distinct elementary disk types with non-zero 
coefficients, four triangles and one quadrilateral. If the j*'' elementary disk type 
occurs, then fix one of the edges that it meets, and add 1 to the j*^ component of 
the weight vector at each of the indices that the j*'' elementary disk meets on that 
edge. The orbit weights are then the normal coordinates of the components of the 
normal surface F. Each point in the output list L' corresponds to a component 
of the normal surface with normal coordinates given by the corresponding weight 
in Z^* . Theorem El tells us that the list L' is computed in time polynomial in 
A:m(ilogDlog A^. We now bound these constants in terms of t. 

Since each edge of a tetrahedron meets at most three disk types in that tetra- 
hedron, each edge of a tetrahedron can contribute at most six points at which the 
weight vector changes. Given six edges to each tetrahedron, we have m < 36t. As 
before we have a bound for the number of pairings k < 12t and the number of 
normal coordinates is given hy d = It. The total weight bounds the normal coor- 
dinates, and with D = W we get a bound on the running time that is polynomial 
intlogW. □ 

Corollary 18. Let M be a 3-manifold with a triangulation T containing t tetra- 
hedra and let F be a fundamental normal surface in M . There is a procedure for 
counting the number of components of F and the topology of each component which 
runs in time polynomial in t. 

Proof. Theorem 01 gives a bound for the normal coordinates of F of t2^*+^. Recall 
that there are It normal coordinates, and each represents a triangle or quadrilateral, 
so the total number of intersections with the 1-skeleton satisfies W < 28t^2^*"''^. In 
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particular, logTy is bounded above by a polynomial in t. Plugging this in for W in 
Corollary 1 1 71 we get a bound for the running time which is a polynomial in t. □ 

7. The complexity of minimal spanning area 

In this section we examine the complexity of the problem of determining the 
smallest area of a spanning surface for a curve in a 3-dimensional manifold. Such 
an area calculation problem seems at first to be ill suited to a complexity analysis, 
since it has real solutions depending on a choice of Riemannian metric. 

We recast the area calculation problem into a discretized form where its complex- 
ity can be analyzed. Given a curve in a suitably discretized Riemannian 3-manifold, 
we ask whether it bounds a surface of area less than C, where C is an integer. To 
describe a metric on a 3-manifold with a finite amount of data, we restrict to piece- 
wise flat metrics, and manifolds constructed from collections of flat tetrahedra and 
triangular prisms whose faces are identified by isometrics. The curvature of such 
PL metrics can be defined as a limit of smooth curvatures, and is concentrated 
along their edges and vertices. A particular manifold in this class is described by 
a decomposition into tetrahedra or triangular prisms with a rational (or integer) 
length assigned to each edge. The metric on this tetrahedron or prism is then taken 
as the metric on the Euclidean tetrahedron or prism with those edge lengths. In 
the case of a prism we also set the angles of quadrilateral faces to be right an- 
gles. Prisms are allowed in this construction in order to form metrics with rational 
lengths on spaces that are products. Identified 2-dimensional faces are required to 
be isometric. We do not require that the total angle around an edge is 27r, nor 
do we make any metric conditions at a vertex. This type of metric is described 
by a finite set of data, and can be used to approximate Riemannian metrics on a 
manifold. Up to scaling, we can take all the edge lengths to be integers. We call 
these objects metrized PL 3-manifolds. A curve is given as a collection of edges in 
the 1-skeleton of M. We will show that given an integer C, determining whether 
the smallest spanning surface for a curve in such a 3-manifold has area less than C 
is NP-hard. 

Problem: MINIMAL-SPANNING-AREA 

INSTANCE: A 3-dimensional metrized PL manifold M, a 1-dimensional 
curve K in the 1-skeleton of Af , and a natural number C. 
QUESTION: Does the curve bound a surface of area A < CI 

The size of an instance is given by the number of bits needed to describe all the 
edge lengths and C. 

Theorem 19. MINIMAL-SPANNINC-AREA is NP-hard. 

Proof. We reduce in polynomial time an instance of the NP-hard problem ONE- 
IN-THREE SAT to an instance of MINIMAL-SPANNING-AREA. This shows that 
MINIMAL-SPANNING-AREA is at least as hard, up to polynomial time reduction, 
as ONE-IN-THREE SAT. 

As a first step, we set up a 2-dimensional version of MINIMAL-SPANNING- 
AREA. We then construct a 3-manifold by a thickening process, with the property 
that a minimizing surface must remain within the 2-complex. 

Given a boolean expression representing an instance of ONE-IN-THREE SAT, 
we construct a triangulated metrized 2-complex and an integer C. This complex 
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contains a curve K with the property that the expression admits a satisfying as- 
signment if and only if K bounds a surface of area less than C. This metrized 
complex is shown in Figure [S] for the expression (xi V X2 V X3) A (xi V 2:2 V X3). 




Figure 5. A metrized branching surface corresponding to the 
boolean expression {xi V X2 V X3) A (xi V a;2 V 2:3). The picture 
is not to scale. The shaded prisms are constructed to each have 
area one, while the rest of the surface has total area less than 1/2. 

The branching surface is similar to the one used in the proof of Theorem ^ but 
carries the additional structure of a metrized triangulation, whose triangles have flat 
metrics of prescribed edge length. The metrized triangles are constructed so that 
near each of the m boundary components corresponding to clauses of the boolean 
expression there are three triangulated disks of area close to one, one on each of 
the three handles coming into the punctured sphere near the boundary component. 
These disks are shaded in Figure El Each of these shaded disks is chosen to have 
area between 1 and 1 + l/2m. The surface is constructed so that the union of all 
triangles in the rest of the surface has total area less than 1/2. 

We saw in Theorem ^ that a spanning surface which has minimal genus goes 
over each of the shaded disks at most once and goes over exactly one shaded disk 
for each of the m clauses. It follows that such a surface has total area m < A < 
TO + 1. Furthermore, a satisfying assignment for ONE-IN-THREE SAT leads to 
an embedded spanning surface, with the satisfying values of the variables selecting 
branches of the surface, and such a spanning surface has area less than m + 1. So 
an instance of ONE-IN-THREE SAT can be reduced to an instance of MINIMAL- 
SPANNING- AREA for this 2-complex. 

To pass to a 3-manifold, thicken each triangle in the branching surface to a 
triangular prism, triangulated as in the proof of Theorem ^ and with a product 
metric. This produces a 3-manifold M which is a thickened up version of the 
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2-coinplex. Projection to the branching surface is area non-increasing, and area 
decreasing for a surface with boundary on the branching surface but not contained 
in it. Therefore a least area surface spanning K must He on the branching surface. A 
closed manifold DM with a piecewise-smooth metric can be obtained by a doubling 
construction as in Theorem ^ The doubling involution is an isometry, so that 
reflecting a surface meeting DM\M into M does not increase area. It follows that 
the embedded spanning surface on the branching surface is a least area surface in 
DM. □ 



8. Open questions 
Among many unresolved questions are: 

1. Does determining knot genus remain NP-hard if we restrict to knots in the 
3- sphere? 

2. Is determining the genus of a knot in a 3-manifold NP? This amounts to 
showing that finding a lower bound to the knot's genus is an NP problem, in 
contrast to the upper bound we have investigated. Recall that the genus of a knot 
is the least possible genus of all spanning surfaces. We have shown that certifying 
that the genus is at most g is NP, but left open the possibility that the genus may 
be smaller than g. If the answer to this question is yes, then we can certify that a 
non-trivial knot has positive genus, and it would follow that UNKNOTTING is 
both NP and coNP. 
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